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Abstract 

We present, in the framework of the interacting hadron resonance gas, an evaluation 
of thermodynamical quantities. The interaction is modelled via a correction for the 
finite size of the hadrons. We investigate the sensitivity of the model calculations 
on the radius of the hadrons, which is a parameter of the model. Our calculations 
for thermodynamical quantities as energy and entropy densities and pressure are 
confronted with predictions using the lattice Quantum Chromodynamics (QCD) 
formalism. 



1 Introduction 



One of the major goals of ultrarelativistic nuclear collision studies is to obtain informa- 
tion on the QCD phase diagram [I]. A promising approach is the investigation of hadron 
production. Hadron yields measured in central heavy ion collisions from AGS up to RHIC 
energies can be described very well [2f3P| 5,6,7,8,9j within a hadro-chemical equilibrium 
(also called hadron resonance gas or statistical) model. The main result of these investi- 
gations is that the extracted temperature values rise rather sharply from low energies on 
towards a center-of-mass energy per colliding nucleon pair ^/snn —10 GeV and reach, for 
higher collision energies, constant values near T= 160- 165 MeV, while the baryochemical 
potential fit, decreases monotonically as a function of energy. The Hagedorn limiting tem- 
perature [TU] behavior suggests a connection to the phase boundary between the hadronic 
world and the deconfined phase. It was, indeed, argued [TT] that the quark-hadron decon- 
finement phase transition drives the equilibration dynamically, at least for SPS energies 
and above. For lower energies, the quarkyonic state of matter [12] could complement this 
picture. The conjecture of the triple point [13] between hadronic, deconfined and quarky- 



onic matter was put forward in this context. In a recent study [H] it is argued, however, 
that the chemical frezee-out region at large /x& is not close to the phase boundary. 

Theoretical investigations of the QCD phase diagram are an important priority for on- 
going research. While effective field theories need to be employed to model QCD in the 
strongly interacting regime for finite /i& [T5"fl6|ll7j . QCD calculations on the lattice are 
an increasingly reliable approach for /i& ~ 0. Employing calculations of QCD on lattice 
all groups predict indeed a steep increase of thermodynamical quantitites near a critical 
temperature for deconfinement, T c . Until recently, values for T c between 151 MeV [18] 
and 192 MeV [19] were obtained. New results on larger lattices and with quark masses 
approaching the physical values [2T)f2"T] lead to better agreement on the T c value in the 
range 155-160 MeV, while important details of lattice QCD calculations continue to be 
addressed. In this context, the hadron resonance gas (HRG) is used by lattice QCD groups 
as reference for their calculations in the hadronic sector [221123] . This follows earlier ideas of 
pM] and involves modelling of quark mass dependence in the hadron resonance gas model 
to account for the finite lattice spacing [2"2l23f24] . Conversely, lattice data are used to 
constrain effective models based on hadronic resonances [23] to describe hadronic matter 
near T c . 

The details of the hadron resonance gas (or statistical) model are important too, as we 
have recently shown in [9], where we demonstrated that the completeness of the hadron 
spectrum involved in calculations is important for a precise description of data in nucleus- 
nucleus collisions. The aim of this paper is to confront our HRG model calculations with 
lattice QCD predictions for thermodynamical observables in the hadronic sector. In par- 
ticular, we investigate the effect of the excluded volume correction employed in the model 
to approximate a short-range repulsive hadron-hadron interaction. Such excluded volume 
corrections were first introduced in [26], albeit not yet in a thermodynamically consistent 
way. A thermodynamically consistent approach was first developed in [27] and will be the 
basis for our investigations. 



2 Model description 



We restrict ourselves here to the basic features and essential results of the statistical model 
approach. A complete survey of the assumptions and results, as well as of the relevant 
references, is available in ref. [2B] . 

The basic quantity required to compute the thermal composition of hadron yields and 
the thermodynamical quantities is the partition function Z(T, V). In the grand canonical 
(GC) ensemble, the partition function for a particle species % in the limit of large volume 
takes the following form (k = h = c = 1): 
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from which the particle density rij, the partial pressure Pi, the energy density S{ and the 



entropy density s, are then calculated according to: 
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where ^ = (2Jj + 1) denotes the spin degeneracy factor, T is the temperature and Ei = 



p 2 + mf is the total energy. The (+) sign corresponds to fermions and (-) corresponds to 
bosons. For the hadron species i of baryon number Bi, third component of the isospin 
strangeness Si, and charm Cj, the chemical potential is = + /// 3 i3j + /isSj + ficCi- 
The chemical potentials related to baryon number (fib), isospin (/ij 3 ), strangeness (fic) an d 
charm (/i^) ensure the (on average) conservation, in the collision, of the respective quan- 
tum numbers: i) isospin: V cons J2i n ihi — ^3°*, with V cons = N B / J2i n iBi] ii) strangeness: 
J^i^iSi = 0; hi) charm: J2i n iCi = 0. The (net) baryon number N B and the total isospin 
I^ ot of the system are input values which need to be specified according to the collid- 
ing nuclei and rapidity interval studied. Taking into account the conservation laws (i-iii), 
the freeze-out temperature T, the baryochemical potential fib and the fireball volume at 
chemical freeze-out V are the parameters of the model, which are obtained from fits to 
experimentally measured hadron yields. 

The following hadrons (number of species, not counting are included in the calcula- 
tions: i) mesons: non-strange (123), strange (32), charm (40), bottom (28); ii) baryons: 
non-strange (48), strange (48), charm (32), bottom (14). The corresponding anti-particles 
are of course also included. Their characteristics, including a rather complete set of decay 
channels (all strong and electromagnetic decays), are implemented according to the 2008 
PDG compilation_j [22], with hadron masses reaching 3 GeV. We use vacuum masses for 
all hadrons. 



Usually, whenever thermal fits are performed, the finite widths of resonances are taken into 
account in the density calculation by an additional integration, over the particle mass, with 
a Breit- Wigner distribution as a weight [8] . For the range of temperatures investigated in 
this work the effect of the finite resonance widths is small and to save computing time 
we have not employed the additional integration in the present calculations except where 
stated otherwise. 



The 2010 PDG compilation contains updates in the hadron mass spectrum, but these are 
expected to have a minor influence on our results. 



3 Interactions in the hadron gas model 



When comparing thermodynamical quantities computed within the framework of the 
hadron resonance gas model with results obtained using lattice QCD methods one has to 
decide how to incorporate interactions among the hadrons. One approach is to use results 
obtained by the authors of [30|31|32j where two-body collisions are taken into account 
through scattering phase shifts. Here the interaction measure (the 2nd virial coefficient) 
is related to the derivative of the phase shifts with respect to energy. To compute the 
thermodynamics of the interacting hadron resonance gas in this way one would need 
knowledge of the energy dependence of all phase shifts. At first glance this seems quite 
impractical. An interesting result was, in this context, obtained in [33j. These authors 
show by explicit construction that, for simple systems such as gases of pions, pions and 
nucleons, and pions, kaons, and nucleons, the equation of state of the interacting system is 
obtained by adding the relevant resonances, the p and /°(980) mesons, the A baryon, the 
K* meson, to the list of particles and by computing the partition function of the enlarged 
gas assuming no interactions. 

This interesting result has led some authors [22 S 2S] to argue that the thermodynamics 
of the interacting hadron resonance gas is well approximated, via the Dashen, Ma and 
Bernstein theorem [3~Tf32] . by that of the non-interacting case, provided that all states 
(resonances) are included in the partition function. It is one of the goals of this paper to 
address the accuracy of this approximation within the framework of our interacting HRG 
model. Even at the formal level, there are points to be considered. First, the t], u, 77', <f> and 
ao mesons cannot be treated like this [33]. Second, the baryon-baryon interaction is largely 
repulsive, with no known resonance structure^] see, e.g. fig. 8 of |33j . More importantly, 
the approach of [3T1I32] is, as also discussed there, a low density approach, relevant for 
dilute systems. At temperatures near T c , the temperature of the phase boundary between 
hadron gas and quark-gluon plasma, the hadron resonance gas is not dilute anymore. As 
will be shown below, overall densities exceed 0.5 fm -3 and total baryon densities are close 
to normal nuclear matter densities of 0.15 fm -3 . This implies that for the whole range 
of baryon chemical potentials considered here the baryon densities near T c are close to 
or exceed the critical value worked out in [33] above which the virial expansion breaks 
down. In this environment also the concept of asymptotic states needed for the S-matrix 
approach of [3"0131II32] is ill defined. 

We therefore explore in the following, in addition to the 'free' hadron resonance gas, also 
the thermodynamic properties of a hadron resonance gas in which short-distance repulsion 
is explicitely taken into account using the thermodynamically consistent excluded volume 
approach developed in [27]. In essence this amounts to a Van-der-Waals construction. This 
is implemented according to [2"T1I3"1] in an iterative procedure for the total pressure as: 



where p td -9 as = Y^iPt ' 9as (T,fti) and for each particle i the chemical potential at a given 
iteration is recalculated as: 
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2 We neglect here the deuteron in the 3 Si state. In the baryon-antibaryon system there is likely 
no short-range repulsion and this leads to a small correction which is discussed below. 



This approach yields the following formulae for the particle densities rij, the total en- 
ergy density e and the total entropy density s expressed as a function of the respective 
quantities in the ideal gas, which are given in (J^JtSJ) ^ 



Hi = rii(T,fii, ...,fi m ) 



dP\ nf 9as (T,/i. 



dnij T 1 + Efe V eigen)k n^ 9as (T, f2 k ) 



s = s(T,^i,...,^ m )= — = — — — ld - (9) 

V 61 ) fll ,..., f i m 1 + Efc V eigen , k n k (T, n k ) 



e = e{T,ti 1 ,...,ti m ) = ; — t ^ ^— - ; ~ (10) 



E^f 9ag (T,/2. 

Et r id. gas i rp 

k Veigen,k n k \-L j /^fc 

where V e igen,i = 4 • 47ri??/3 is the eigenvolume of a hadrorf 3 ! with radius We checked 
numerically that thermodynamical consistency, expressed by e = Ts — P + J2i ^i n ii is wen 
fulfilled by our calculations, explicitly confirming the consistency of the procedure 
used for the excluded volume correction. 



For the radius parameter Ri, governing the excluded volume calculation, we follow the 
earlier arguments in [I]. There it was argued that, for baryons, the radius is given by the 
hard-core repulsive interaction as extracted from nucleon-nucleon scattering [3S], giving 
a radius of about 0.3 fm. Values for other baryons should be similar. For mesons, in the 
absence of detailed information on their interactions at short distance, we assign the same 
radius value, based on the similarity of the meson charge radii compared to baryons and 
on the energy dependence of the pion-nucleon phase shifts [36] • For illustration, we have 
included the case of R=0 for mesons, although we believe that the physical case is for 
mesons with repulsive core radius R comparable to baryons. 

For the baryon-antibaryon system there are likely no short-range repulsive interactions, 
because of annihilation processes, which are by construction included in the hadron res- 
onance gas at equilibrium. We have modelled the absence of short-range repulsion in 
a schematic way, introducing a correction factor k ann i for the (anti)baryon eigenvolume 
based on the expression: 

1 _ 1 2Tij 3ar y 0nS 7l an ^j)-| ar y 0nS 

^anni - 1 ~ J— — ~ ^, 

V ' b baryons ~ 1 '•anti baryons J 

where ^baryons is the density of baryons and n ant ibaryons the density of antibaryons. As 
shown by the results presented in the following, the absence of short-range repulsion in 
the baryon-antibaryon system leads to only a small correction, since mesons dominate at 
small fif, and since there are very few antibaryons at large 

For the rest of the paper we show results of our calculations for radii in the range of 
0.3±0.05 fm, along with the cases discussed above and contrast those with the case of no 
interaction (R4—O). The value of 0.3 fm, common for mesons and baryons, was also used 
whenever we performed thermal fits to hadron abundancies [4~f6|8|l9"] . In the description of 
hadron yields with the statistical model, the excluded volume correction leads to a larger 
volume parameter, while the fit temperature and baryochemical potential are unchanged 

3 Consider a particle with radius R in the hard sphere model. Then no other particle can 
come closer than a distance 2R. Per pair the excluded volume is 47r(2i?) 3 /3, leading to V e i gen = 

4 • 47ri? 3 /3 for the particle. 



(11) 




Figure 1. Total hadron densitites as a function of temperature, calculated for the hadron reso- 
nance gas model at /i&=0, without (dashed line) and with (band for R=0.3±0.05 fm) the excluded 
volume correction. The case of -R m eson=0 is shown with the dotted line, while the dot-dashed 
line denotes the effect of the baryon-antibaryon annihilation correction. 

compared to the case of fitting hadron ratios [8], for the case of identical Ri for mesons 
and baryons. The implication of a pion radius different from all other hadrons for the 
description of data has been studied by Yen et al. [34] . 

In Fig. [1] we present the temperature dependence of the hadron densities calculated with 
our model, with and without interactions modelled via excluded-volume corrections, as 
well as with the cases of i? m cson=0 and of absence of short-range interactions for baryon- 
antibaryon pairs. This illustrates our remarks above, namely that, while at low tempera- 
tures (low densities) there is no difference between the case of excluded volume correction 
and the case of free hadron resonance gas, the difference becomes appreciable for T above 
130-140 MeV, which is significantly below the value for the critical (crossover) temper- 
ature T c ~160 MeV. Near T c , the hadron resonance gas becomes manifestly dense, with 
the mean distance between hadrons getting significantly smaller than twice the hadron 
radius. All approximations appropriate for the dilute gas, discussed above, break down 
and the non-interacting hadron gas is not a suitable approach anymore. 

We illustrate, for the energy density, e, in Fig. [2] the sensitivity of the calculations on 
the radii of the eigenvolume for mesons and baryons, -R meS on and -Rbaryon, respectively. 
The calculations have been performed for T= 164 MeV, corresponding to the limiting 
temperature reached in heavy-ion collisions [9] and for /i&=0.8 MeV, the value expected 
for the LHC energy according to ref. [9]; calculations for /ife=0 lead to identical results. 
We observe a strong influence of the excluded volume correction on the energy density 
and this is the case for all other thermodynamical quantities. Due to the larger abundance 
of mesons (and in particular of pions) in the hadron gas at these values of T and /if,, the 
sensitivity on R ma thrmmeson is more pronounced. 



Figure 2. The energy density of a hadron gas as a function of the radius (for the eigenvolume 
calculation) for mesons and baryons, at a temperature value of 164 MeV and //{,=(). 8 MeV. The 
dot indicates the radius value of 0.3 fm which we employ as default. 




Figure 3. The adiabatic speed of sound (in units of velocity of light, c=l) squared as a func- 
tion of temperature, calculated for the hadron gas model, with and without excluded-volume 
corrections. Our calculations are compared to lattice QCD calculations of Borsanyi et al. [37] . 



As pointed out earlier [27|33] , a possible problem of the hadron gas model with excluded 
volume corrections is acausal behavior (speed of sound larger than velocity of light). As 
we demonstrate in Fig. [3X our model is not plagued by such a behavior (for the case /15 ~0 



considered here). The adiabatic speed of sound, c s , is calculated as: 




dlnT 



din s 
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It exhibits a shallow minimum as a function of T for the case of excluded volume correc- 
tions. Our calculations are compared to lattice QCD calculations [37] • In general, a good 
agreement between our and the lattice result is observed. We note that our calculations 
predict a shallow minimum for T around 140-150 MeV, while the lattice values indi- 
cate a more pronounced dip and exhibit a speed of sound value smaller than the hadron 
resonance gas with interactions. It would be interesting to see if the corresponding low 
temperature part of the equation of state (EoS) would lead to changes in hydrodynamic 
calculations, where generally the EoS shown in Fig. 6 or ref. [22] is used. 



4 Hadron resonance gas and lattice QCD results 

In the following we compute, in the HRG model, thermodynamical quantities with and 
without excluded volume corrections and compare the results to predictions from lattice 
QCD. In Fig. H] we show the temperature dependence of energy density, pressure, and 
entropy density, each normalized to appropriate powers of the temperature. 

The case without interactions (no excluded volume correction) has the expected strong 
dependence on temperature. As noted early on by Hagedorn [ID] , a limiting temperature, 
also called "Hagedorn temperature", of Tjj ~ 200 MeV arises for calculations of ther- 
modynamical quantities within the HRG model if one assumes a hadron mass spectrum 
which increases exponentially with particle mass. Such exponential behavior is consistent 
with the present knowledge of hadron resonances [29] up to 2.0-2.5 GeV in mass. At 
this temperature all thermodynamical quantities for the HRG without excluded volume 
corrections diverge. We note in passing that all thermal model calculations without ex- 
cluded volume corrections become meaningless for temperature values close to Tjj. In the 
course of investigations reported in [9] we realized that this implies a practical limitation 
to temperatures below 175 MeV as all calculations for higher temperatures become very 
sensitive to details of the mass spectrum for masses larger than 3 GeV. 

For the case of calculations employing finite hadron volume corrections the Hagedorn 
infinities are tamed. This was already noted by Hagedorn [10] who was the first to intro- 
duce excluded volume corrections [25]. Our findings substantiate this and imply that the 
Hagedorn limiting temperature is an artifact of the usage of the free hadron resonance 
gas description at temperatures where the implicit approximations for dilute systems are 
manifestly inappropriate. 

For temperatures below 120 MeV the HRG model results with and without excluded 
volume correction almost coincide, see Fig. HI For larger temperatures, the HRG with 
interactions yields, in our view, a realistic description of the hadronic phase. Therefore, in 
the confined regime, lattice QCD calculations of thermodynamical variables should give 
results in agreement with the interacting HRG. The expectation is that, as soon as effects 
of deconfinement become important in the lattice QCD results, they should increasingly 
exceed the HRG values. In Fig. H]the most recent predictions of lattice QCD are compared 




Figure 4. Temperature dependence of thermodynamical quantities. The calculations with the 
hadron gas model are shown without (dashed line) and with (band for R=0.3±0.05 fm) the 
excluded volume correction. The case of i? m eson=0 is shown with the dotted line, while the 
dot-dashed line denotes the effect of the baryon-antibaryon annihilation correction. They are 
compared to LQCD results of Borsanyi et al. [37] . 

to the HRG results. Indeed, below T=150 MeV good agreement between results of lattice 
QCD [37] and the interacting HRG is found. On the other hand, effects of the onset of 
deconfinement [201I2T] are apparent for T in excess of 150 MeV. 

In Fig. [5] we compare our calculations for the trace anomaly e — 3P (normalized to T 4 ) 
to LQCD results [37f38] . We see an agreement between LQCD data and our calculations 
for the interacting hadron gas only up to T=140 MeV for the data of Borsanyi et al. [37] . 
For the (preliminary) data of the HotQCD collaboration [35] we have used for illustration 
the set for hisq action with N t =Q, but we note that other available sets are in agreement 
with those within the errors [38] . 

We turn now to the energy dependence of the thermodynamical quantities at chemical 
freeze-out in central nucleus- nucleus collisions [9]. The degree of stopping of the colliding 
nuclei, which is energy dependent, brings some uncertainty in the choice of Nb and I^ ot . As 
we study central collisions of heavy nuclei (Au or Pb) and focus on data at midrapidity, 
we have chosen Nb = 400 • /i 6 /938 MeV and If* = -40 ■ /i fe /938 MeV, reflecting that 
fib traces stopping of the two colliding nuclei. The sensitivity of the thermodynamical 
quantities on Nb and I^ ot is, however, small. 




T (MeV) 



Figure 5. Temperature dependence of the trace anomaly. The calculations within the hadron 
gas model (lines, as in Fig. H|) are compared to LQCD calculations of Borsanyi et al. [37] and 
HotQCD collaboration [38] (preliminary results). 

In Fig. [6] we show the thermodynamical quantities as a function of collision energy for 
chemical freeze-out in central nucleus-nucleus collisions. The trends seen in Fig. [6] re- 
flect primarily the sharp increase of the temperature at chemical freeze-out (determined 
from fits of experimental data up to v /J/v r /v=200 GeV [9]) followed by a saturation above 
v /J/v r /v —10 GeV. The characteristic energy dependence of the baryon density, exhibiting a 
maximum around 8 GeV, is determined by the increase of T combined with the strong de- 
crease of fib with energy, as discussed in [9] . The effect of interactions leads to up to 30% 
reduction of the thermodynamical quantities at chemical freeze-out in nucleus-nucleus 
collisions. 



5 Summary 

We have presented an evaluation of thermodynamical quantities in the framework of the 
interacting hadron gas model, incorporating all known hadrons with masses reaching 3 
GeV. A Van-der-Waals-type interaction is modelled via an excluded volume correction. 
Thermodynamic consistency is ensured by construction and the model exhibits proper 
causal behavior. 

The resulting values for the thermodynamical quantities increase, already for temperatures 
of 130-150 MeV, i.e. significantly below T c , much less steeply than in case of a free hadron 
resonance gas. Near T c the free hadron resonance gas calculations already show signs 
of the Hagedorn divergence. Comparisons of lattice QCD results with the free hadron 
resonance gas in this temperature regime are therefore in our view problematic. Our 




Figure 6. Energy dependence of energy density, pressure, entropy density and baryon and meson 
densities at chemical freeze-out in central nucleus-nucleus collisions. The full lines are for the 
excluded volume corrections, the dashed line without, while the case of R m athrmmeson = ^ is 
shown with the dotted line and the dot-dashed line denotes the effect of the baryon-antibaryon 
annihilation correction. 



results imply the need to consider the hadron resonance gas with interactions, beyond the 
usual implementations based on the Dashen, Ma and Bernstein theorem. 

On the other hand, lattice QCD simulations start to be precise enough to reproduce the 
complete hadron gas at low temperatures. The apparent rise of the lattice QCD results 
above the HRG results is a clear indication of the onset of deconfinement not contained 
in the latter. In our view, the lattice results show genuine quark and gluon degrees of 
freedom in the vicinity of the (crossover) transition. In this temperature range, the lattice 
results produce thermodynamical quantities well above our predictions for the interacting 
hadron resonance gas. Our findings also imply that the Hagedorn limiting temperature 
is an artifact of the usage of the free hadron resonance gas description at temperatures 
where the gas becomes manifestly dense. 
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